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Abstract. The helium-atmosphere (DB) white dwarfs are commonly thought to be the descendants of the hotter PG 1159 
stars, which initially have uniform He/C/O atmospheres. In this evolutionary scenario, diffusion builds a pure He surface layer 
which gradually thickens as the star cools. In the temperature range of the pulsating DB white dwarfs (T c g ~ 25 000 K) 
this transformation is still taking place, allowing asteroseismic tests of the theory. We have obtained dual-site observations of 
the pulsating DB star CBS 114, to complement existing observations of the slightly cooler star GD 358. We recover the 7 
independent pulsation modes that were previously known, and we discover 4 new ones to provide additional constraints on the 
models. We perform objective global fitting of our updated double-layered envelope models to both sets of observations, leading 
to determinations of the envelope masses and pure He surface layers that qualitatively agree with the expectations of diffusion 
theory. These results provide new asteroseismic evidence supporting one of the central assumptions of spectral evolution theory, 
linking the DB white dwarfs to PG 1 159 stars. 
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1. Introduction 

The origin and evolution of the helium-atmosphere (DB) white 
dwarf stars has long been a subject of debate. They com- 
prise roughly 20% of the population of field white dwarfs, 
with most of the remaining 80% consisting of their hydrogen- 
atmosphere (DA) cousins. The majority of white dwarfs in both 
of these spectral classes are thought to arise from the evolu- 
tion of isolated main-sequence stars with masses that are in- 
sufficient to ignite carbon fusion, which ultimately leave their 
hot carbon/oxygen cores to descend the white dwarf cooling 
nack. The bifurcation into two spectral classes is thought to 
occur during post-asymptotic-giant-branch (post-AGB) evolu- 
tion when, in some cases, a very late thermal pulse burns off the 
residual hydrogen in the enve lope, producing a nearly pure he- 
lium atmosphere dlben et al.l l 1983 >. Such objects are then sup- 
posed to return to the white dwarf cooling track as PG 1159 
stars, which are widely believed to be the precursors of most 
DB white dwarfs. 

The primary difficulty with this scenario is the paucity 
of DB white dwarfs with temperatures between 45 000 
and 30 000 K — a phenomenon known as the "DB gap" 
dLiebert et al.lll986l) . which persists even in the latest re sults 
from the Sloan Digital Sky Survey dKleinman et al.ll2004 . If 
the PG 1159 stars and DB white dwarfs are linked, then why 
don't we see the long-lived intermediate phase? The spec- 



tral ev olution theory advanced bv lFontaine & Wesemaelldl987l 
Il997l) postulates that traces of hydrogen (M H ~ 10" 15 M«) in 
the envelopes of hot PG 1159 stars float to the surface and 
disguise them as DA stars within this temperature interval. 
A growing helium convection zone eventually dilutes the hy- 
drogen in the photosphere until it is no longer detectable at 
30000 K, and the star continues its evolution as a DB. A rel- 
ative overabundanc e of DA stars inside th e DB gap supports 
this hypothesis (see Kleinm an et alJl20041 their Fig. 11). But 
lower than predicted surface hydrogen abundances, measured 
for several DB stars just below the gap, remain a problem for 
spectral evolution theory dProvencal et al J2000t> . 

If we assume that there is an evolutionary connec- 
tion between the PG 1159 stars and the cooler DB white 
dwarfs, we can look to several independent groups who 
have used time-dependent diffusion calculations to follow the 
chang es in the interior structu r e of these objects as they 
cool dDehner & Kawaletl 1 19951: IFontaine & BrassardI |2002; 
lAlthaus & CorsicolkOO^ The hot PG 1159 stars, having re- 
cently emerged from the born-again phase, contain envelopes 
with a nearly uniform mixture of heliu m (He), carbon (C), and 
oxygen (O) out to the photosphere dDreizler & Heberl fl998: 
Her wig et aljfl 999). As they cool, the He diffuses upward and 
gradually accumulates to form a chemically pure surface layer. 
This leads to a double-layered structure, with the pure He sur- 
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face layer overlying the remainder of the uniform He/C/O en- 
velope, all above the degenerate C/O core. 

A key prediction of the diffusion models is that, for a 
given stellar mass, the pure He surface layer will steadily 
grow thicker as the DB star cools. The only available observa- 
tional tests of this prediction come from asteroseismology — the 
study of the internal structure of stars through the interpreta- 
tion of their pulsation periods. Fortunately, non-radial oscilla- 
tions are naturally excited in otherwise normal DB white dwarf 
stars as they cool through a narrow range of temperatures near 
25 000 K. However, there are only 9 such pulsating DB stars 
(DBVs) presently kn own, wi th effective temperatures between 
22 400 and 28 400 K jBeauchamp et al.ll999l) . 

The most thoroughly studied DBV star is GD 358, which 
has been the target of three coordinate d observing cam- 
paigns of the Whole Earth Telescope (WET:lNather et alJl990t 
IWinget et alJ 1 1994 IVuille et alJ 1200 0: Kepl er et all 120031) . 
These observations revealed dozens of independent pulsa- 
tion periods, inclu ding a sequence of 11 dipole (6=1, m=0; 
KoIak_gLaD I2QQ3J) mode s of consecutive radial overtone (k; 
Bradl ey & Wing ed [l994), each providing a complementary 
probe of the stellar interior. These data formed the basis for 

the only asteroseismic tests of the diffusion models to date 
1 r II I 1 II I. 

(Dehner& Kawaler 1995; Fontaine & Brassard 2002). In the 

most extensive of these tests, Fontaine & Brassard (2002) pro- 
duced a targeted grid of double-layered evolutionary models to 
try to fit the periods of GD 358 in detail. They found evidence 
of two composition transition zones near outer mass fractions 
of 10~ 5 8 and 10~ 3 M„, which they interpreted as the base of 
the pure He surface layer and th e still-uniform He/C/O enve- 
lope, respectively. However, Montgom ery et alJ 1120031) showed 
that there is a potential ambiguity in the inferred locations of 
these internal structures, due to an inherent core-envelope sym- 
metry for models of DBV stars. As a consequence, real struc- 
ture at specific locations in the core could be mis-interpreted as 
structure at symmetric points in the envelope, and vice versa. 

The best way of circumventing this potential ambiguity 
is to fit additional stars using the same models, to determine 
whether the independent observations are consistent with the 
same underlying physical processes. Such model-fitting at- 
tempts for DBV stars other than GD 358 have been hampered 
by the lack of suitable observational data. Handler et al. (2002, 
hereafter HMW) found a total of 7 independent puls ation peri- 
ods in the slightly hotter DBV star CBS 1 14, and Met calfe et alJ 
(2003) attempted to fit these observations using double-layered 
envelope models that were also applied to GD 358. They found 
qualitative agreement with the predictions of diffusion theory 
(thinner surface He layers for hotter white dwarfs of compara- 
ble mass) when using models that also included an adjustable 
C/O profile in the core. However, the large number of free pa- 
rameters in their model, relative to the number of observed pe- 
riods in CBS 1 14, made these results somewhat uncertain. 

2. Observations & Data Reduction 

Motivated by the need for additional pulsation periods to con- 
strain the model-fitting of CBS 1 14, we organized a dual-site 
campaign using 2 m class telescopes in February 2004. Our 



Table 1. Journal of observations. 



Run 


Telescope 


Date 


Start 


Length 








(UT) 


(s) 


040219 


BOAO 1.8m 


2004 Feb 19 


11:15:50 


33540 


040220 


BOAO 1.8m 


2004 Feb 20 


10:58:19 


13080 


A0844 


McD 2.1m 


2004 Feb 21 


04:17:05 


11050 


040223 


BOAO 1.8m 


2004 Feb 23 


10:42:57 


29520 


040224 


BOAO 1.8m 


2004 Feb 24 


15:58:30 


14520 


040225 


BOAO 1.8m 


2004 Feb 25 


12:49:33 


24840 


A0847 


McD 2.1m 


2004 Feb 26 


01:53:07 


36120 


A0848 


McD 2.1m 


2004 Feb 27 


01:51:17 


38550 



primary goals were to confirm the periods of the 7 pulsation 
modes found by HMW, and to discover additional modes below 
the ~5 milli-magnitude noise limit of the previous observations 
(obtained over 3 weeks from the 0.75 m telescope at SAAO). 
This campaign was awarded 7 nights (2004 Feb 19-25) on the 
1.8 m telescope at Bohyunsan Optical Astronomy Observatory 
(BOAO) in Korea, and 7 nights (2004 Feb 21-27) on the 2. 1 m 
telescope at McDonald Observatory (McD) in Texas. 

At BOAO we used a thinned SITe 2K CCD camera with 
no filter. To minimize the readout time, we binned the pixels 
2x2 for an effective resolution of 0.68 arcseconds pixeF 1 , and 
we used only the central portion of the detector for a field of 
view ~4 arcminutes on a side. On the first night we used an 
exposure time of 8 seconds for a cycle time of 20 seconds, but 
on subsequent nights we increased the exposure time to 20.5 
seconds for a cycle time of 30 seconds. This system has a gain 
of 1.8 <? ADU 1 and a readout noise of 7 e~ rms. 

The observations at McD were obtained using the Argos 
time-series CCD photometer dNather & Mukadam 2004) with 
no filter. This instrument uses a back-illuminated 512x512 
frame-transfer CCD mounted at the f/3.9 prime focus, for a 
field of view 2.8 arcminutes on a side and a resolution of 0.33 
arcseconds pixeF 1 . We used an exposure time of 10 seconds, 
which was also the cycle time — the transfer of each image to 
the masked on-chip buffer requires only 3 10 /is, and the 280 ms 
readout takes place during the next exposure. This system has a 
gain of 2 eADV 1 and a readout noise of 8 e~ rms. The obser- 
vations collected for this dual-site campaign are summarized in 
Tabled 

The raw CCD frames were analyzed on-the-fly during the 
observing campaign to monitor the data quality, but final reduc- 
tions were done uniformly using the IRAF high-speed photom- 
etry (HSP) package developed by Antonio Kanaan (personal 
communication). After correcting each image for bias level and 
flat field effects, we extracted aperture photometry for CBS 1 14 
and three nearby comparison stars. The aperture was typically 
set to have a radius between 2 and 4 arcseconds with a sur- 
rounding annulus to determine the local sky level, but the aper- 
tures were adjusted for each data set to maximize the S/N of 
the resulting light curves. The time-ser ies measurements were 
then imported into XQED (Riddle 2003), the photometry anal- 
ysis software for the WET. For each light curve, we subtracted 
the local sky contribution and performed a point-by-point di- 
vision of the target star by the average of the three compari- 
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Fig. 1. Reduced light curves for CBS 1 14 during our dual-site campaign in 2004, showing fractional intensity variations of +15 
percent with a dominant period of ~650 seconds. 



son stars. We then corrected for any low-frequency variations 
(e.g., caused by differential color extinction) by fitting and re- 
moving a second-order polynomial. Finally, we subtracted the 
mean light level and applied barycentric timing corrections to 
the exposure mid-times. The resulting light curves for CBS 1 14 
are shown in Fig.^ 

3. Frequency Analysis 

The Fourier Transforms (FTs) of these light curves, for each 
observatory and for the combined data set, are shown in Fig. 13 
The corresponding window functions — the FT of a single si- 
nusoid sampled in the same way as the data — are shown to 
the same scale in each panel. The distributed sampling of the 
BOAO data yields better effective frequency resolution, while 
the slightly larger aperture and higher cadence of the McD data 
leads to a higher S/N. The full data set combines all of these 
strengths, and reduces the amplitude of the 1 d aliases (sepa- 
rated from the main frequency by +1 1.6 //Hz). 

To determine the significance l evel of peaks i n each of these 
FTs, we followed the analysis of Kepler ( 1993). Since power 



is the square of the amplitude, the average power out to the 
Nyquist frequency is a measure of the mean-square deviation 
from zero amplitude, which is dominated by the noise since the 
signal is confined to a relatively narrow range of frequencies. 
The square-root of this value is the root-mean-square ampli- 
tude of the noise, which is analogous to the standard deviation 
(o"noise)- By adopting a significance cut at 3o" no i se , we ensure 
that only ~0.3% of the peaks with this amplitude will not be 
real signals (assuming the noise is Gaussian). Peaks with am- 
plitudes much above the significance cut, and those detected in 
multiple independent data sets, have an even lower chance of 
being due to noise. The 3cr no j se significance level is indicated by 
the dotted line in each panel of Fig. [2] which is at 3.52, 2.33, 
and 2.26 mma for the BOAO, McD, and combined data sets re- 
spectively. For the range of frequencies shown, we expect only 
~4 noise peaks to appear above these cuts. 

We performed a frequency analy sis of the combined data 
set using the program PERIOD98 (Sperl 1998). We started 
by iteratively fitting and subtracting the frequencies with the 
largest amplitudes, refining our estimates of the period, am- 
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Fig. 2. Fourier Transforms (FTs) of the light curves shown in Fig.n° ver the range of independent frequencies exhibited by 
CBS 1 14. The three panels show the FTs for the data from BOAO (top), McD (middle), and for the combined data set (bottom), 
along with the corresponding window functions (inset) and significance cuts (dotted line). 



plitude, and phase for each signal using simultaneous multi- 
ple sine-wave fitting. This allowed us to recover the 7 inde- 
pendent pulsation frequencies found by HMW, with a few im- 
portant differences. We resolved both of the two largest ampli- 
tude modes near 1520 and 1610 //Hz into triplets, though the 
splitting is close to the 1 d aliases (perhaps explaining why the 
earlier single-site data failed to detect them, despite a longer 
time baseline). We also detected an additional component in 
the 1520 //Hz multiplet at 1514.1 //Hz, which could be caused 
by amplitude and/or frequency modulation during the observa- 
tions or might be caused by higher-order combination frequen- 
cies, comparable to the "k = 15" mode seen in GD 358 in 1994 
JVuille et al. 2000). Our frequencies for the modes near 1830 
and 1960 //Hz are both systematically lower than in HMW, but 
by less than would be expected from 1 d alias problems. This 
might be the result of an underlying multiplet structure with 
intrinsic amplitude modulation among the various components 
(e.g., see iRobinson et alJll976l their Fig. 2). We also found 
that our mode near 2300 //Hz was the +1 d alias of the fre- 
quency found by HMW — a possibility that they could not rule 
out. Finally, we resolved the mode near 2510 //Hz, as well as 



a newly discovered mode near 2090 //Hz, into doublets with a 
splitting comparable to the frequency shifts seen in the 1830 
and 1960 //Hz modes. 

In addition to recovering the 7 pulsation modes found by 
HMW, we also discovered 4 new ones. The largest is the dou- 
blet near 2090 //Hz, which nicely fills the gap in the sequence 
of consecutive radial overtones identified by HMW. A marginal 
pea k near this frequency also appeared in the discovery data of 
Winget & Craveini988[ll989» . but it was below the detection 
threshold of HMW. The three other new modes appear at fre- 
quencies near 1190, 1250, and 1380 //Hz. The period spacing 
of these peaks suggests that they are additional independent 
pulsation modes with higher radial overtones than those pre- 
viously seen. Two additional frequencies had amplitudes near 
our significance threshold. These could tentatively be identified 
as an even higher radial overtone near 1 140 //Hz, and a possi- 
ble rotationally split component of the mode near 1720 //Hz. 
These marginal detections prom pted us to com pute weighted 
FTs, following the suggestion of iHandlerl J2003I) . The results 
were nearly identical to the unweighted case, and the marginal 
nature of these two frequencies was unchanged. 
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Fig. 3. The Fourier Transform of our combined data set for CBS 1 14 before and after removing the 28 sinusoids listed in Tabled 
The independent and linear combination frequencies are indicated by the solid and open triangles respectively. The frequencies 
of our optimal model matched to the observed modes are shown as solid squares, while open squares denote model frequencies 
with no detected counterpart. 



Aside from the independent pulsation modes observed in 
CBS 1 14, we also detected 1 1 significant linear combination 
frequencies. Fundamentally, these frequencies arise in our anal- 
ysis because the Fourier Transform decomposes the signal into 
sinusoidal components. A casual inspection of the light curves 
shown in Fig. [^reveals that they are not perfectly sinusoidal. 
Instead, they tend to have sharp peaks and shallow troughs, 
which is commonly interpreted to be associated with a non- 
linear respons e of the surface convection zone to the pulsa- 
tions (e.g., seelBrickhil]|l99llWul200llllsing & KoesterkOOll 
lMontgomervll2005l) . The earlier observations by HMW re- 
vealed two such frequencies, both of which are among those 
we detected. The combination frequencies identified in our data 
involve the six modes or multiplets with the largest amplitudes. 

Our final simultaneous fit of 28 sinusoids to the data using 
PERIOD98 contains a total of 17 independent modes and 11 
linear combination frequencies, which are all listed in Table|2] 
In Fig. |3 we show the FT of our combined data set before and 
after the removal of these 28 frequencies (indicated by solid 
and open triangles for the independent and combination fre- 
quencies, respectively). The marginal signals near 1 140, 1514, 
and 1720 fiHz are apparent in the residuals. 

4. Model Fitting 

To investigate whether these new observations can fit com- 
fortably with the predictions of diffusion theor y, we repea ted 
the glo bal model-fitting procedure described by Met calfe et alJ 
(2003). We used an updated version of the code that incor- 
porates the OPAL radiative op acities d lglesias & Rpgersfl 996) 
rather than the older LAO data ( Huebn er et al.l 19771) . which are 
known to produce systematic errors in the derived temperatures 



Table 2. The multifrequency solution for CBS 1 14. Also shown 
is the mode identification from our optimal model. 

Sinusoid Frequency Period (Amp) Identification 





0"Hz) 


(s) 


(mmag) 


k 


£ m 


Si 


1189.9 


840.38 


3.1 


20 


1 


h 


1252.4 


798.46 


4.3 


19 


1 


h 


1383.0 


723.06 


3.0 


17 


1 


/. 


1509.0 


662.67 


7.5 


15 


1 -1 


St 


1518.8 


658.44 


11.8 


15 


1 


A 


1530.7 


653.30 


9.0 


15 


1 +1 


fi 


1601.2 


624.55 


4.4 


14 


1 -1 


A 


1612.1 


620.31 


15.4 


14 


1 


A 


1622.7 


616.27 


13.3 


14 


1 +1 


/to 


1729.8 


578.09 


7.0 


13 


1 


fu 


1829.6 


546.56 


3.0 


12 


1 


fn 


1963.2 


509.37 


5.5 


11 


1 


Jr. 


2086.5 


479.27 


4.8 


10 


1 ? 


/l4 


2093.6 


477.65 


4.9 


10 




fis 


2317.7 


431.45 


8.9 


9 


1 


/l6 


2510.0 


398.41 


6.8 


8 




/l7 


2515.2 


397.59 


3.2 


8 




/5+/ 6 ... 


3049.4 


327.93 


2.2 






/5+/ 8 ... 


3130.9 


319.40 


3.7 






/5+/ 9 ... 


3141.4 


318.33 


3.5 






ft + ft ■ ■ ■ 


3224.2 


310.15 


2.4 






Ss + fio ■ ■ 


3248.6 


307.83 


2.2 






ft + SlO ■ ■ 


3341.9 


299.23 


3.2 






A + /l2 ■ ■ 


3585.9 


278.87 


2.9 






/5+/15-. 


3836.5 


260.65 


2.4 






f9 + SlS ■ ■ 


3940.4 


253.78 


3.3 






Ss + /l6 • ■ 


4028.8 


248.22 


2.0 






Ss + /l6 • ■ • 


4122.1 


242.59 


2.1 
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Table 3. Optimal model parameters for CBS 1 14 and GD 358. 



Parameter 


CBS 114 


GD358 


Error 


T e ff (K) 


25 800 


23 100 


+ 100 


M,(M ) 


0.630 


0.630 


±0.005 


log(M cnv /M,)... 


-2.42 


-2.92 


+0.02 


log(M Hc /M,)... 


-5.96 


-5.90 


±0.02 


Op (s) 


2.33 


2.26 





jFontaine & Brassard! 19941). The fitting procedure use s a paral- 
lel genetic algorithm ( Metca lfe" & Charbonneaul2003l) to mini- 
mize the root-mean-square residuals between the observed and 
calculated periods (cr P ) for models with effective temperatures 
(Jeff) between 20 000 and 30 000 K, and stellar masses (M«) be- 
tween 0.45 and 0.95 M Q . It allows the base of the uniform He/C 
envelope to be located at an outer mass fraction log(M env /M,) 
between -2.0 and -4.0. The base of the pure He surface layer 
can assume values of log(MH e /A^*) between -5.0 and -7.0. 
To ensure a smooth transition between the self-consistent cores 
and the static envelopes in our models, we fixed the core com- 
position to be pure C. We cannot presently include oxygen in 
the envelopes of our models. 

We applied this fitting procedure to our new observations 
of C BS 114 and to the WET observations of GD 358 from 
1990 dWinget et all 1994 for comparison. For CBS 1 14 we as- 
sum ed that each sing le mode had an azimuthal order m = 
(see Metcalfe 2003a), and we used the central frequency for 
each of the observed triplets. For the two doublets, we used 
the average of the two observed frequencies. The resulting sets 
of optimal model parameters for these two stars are shown in 
Table with internal errors for each parameter set by the res- 
olution of our search. The theoretical pulsation frequencies for 
CBS 1 14 are shown in Fig.[3]as solid squares where they match 
an observed frequency, and the corresponding mode identifica- 
tions are listed in Table |2] For reference, the open squares in 
Fig. show additional model frequencies with no clearly de- 
tected counterpart. 

The spectro s copic temperature determinations of 
iBeauchamp et all i 19991) suggest that CBS 114 could be 
as much as 1500 K hotter than GD 358. Our model-fitting re- 
sults suggest a larger temperature difference, placing CBS 1 14 
at 2700 + 200 K hotter than GD 358. On an absolute scale, 
the temperature and mass of our fit to CBS 1 14 are consistent 
with the spectroscopic values, assuming that any traces of H in 
the atmosphere are small. Our fit to GD 358 is consistent with 
the spectroscopic mass, but our temperature determination is 
slightly low. Our value is also low compared to a more recent 
and precise determination using UV spect roscopy from the 
Hubb le Space Telescope (24 100 + 400 K; Castan heira et alJ 
2005). 

Rec ent time- d epend ent diffusion calculations by 
lAlthaus & Corsicol J2004 suggest that for a given stellar 
mass and envelope thickness, the pure He surface layer should 
thicken by -0.15 dex between the spectroscopic temperatures 
of CBS 1 14 and GD 358. We find a marginal thickening of only 
0.06 + 0.04 dex. However, diffusion models with more massive 
envelopes lead to thicker surface He layers more quickly as the 



star cools, since there is a larger reservoir of He to draw from. 
Our fit to CBS 1 14 has a slightly larger total envelope mass, 
which would tend to diminish the expected difference between 
the thickness of the surface He layers for the two stars to ~0. 12 
dex. The limited grid of published diffusion models currently 
makes more quantitative comparisons impossible. 

5. Discussion 

Our new dual-site observations of the pulsating DB white dwarf 
CBS 1 14 have finally revealed enough independent modes to 
bring the model-fitting for this star to the same level of reliabil- 
ity as for GD 358. We used our updated double-layered enve- 
lope models to match the new observations of CBS 1 14, and the 
archival observations of GD 358 — fitting 1 1 pulsation modes 
in each star. The resulting pair of optimal models have pure He 
surface layers and total envelope masses that qualitatively agree 
with our expectations from diffusion theory. This provides new 
asteroseismic evidence supporting one of the central assump- 
tions of spectral evolution theory, linking the DB white dwarfs 
toPG 1159 stars. 

In addition to the new observational constraints for 
CBS 114, we improved our analysis by incorporating the 
OPAL radiative opacities into the models. Our results include 
high er temperatures and lower masses than the previous fits 
bv lMetcalfe et alJ d2003l) . leading to better agreement with the 
spectroscopic estimates for both stars. The optimal locations 
for the base of the pure He surface layer and the uniform He/C 
envelope are systematically deeper for GD 358 relative to the 
previous fit. For CBS 1 14 the fit is now consistent with the en- 
velope mass e xpected from simulations of carbon dredge-up i n 
DQ stars (see Pell etier et al.ll986tlFontaine & Brassardl2002h . 
and with the He lay er mass e xpected f rom tim e-dependent dif- 
fusion calculations ( Althaus & Corsico 2004). 

These result s have now met the challenge posed by 
Metcalfe (2003b), to fit the observations of more than one DB V 
star and demonstrate that they are both consistent with the same 
underlying physical process. In that paper, the relevant process 
was the nuclear burning history during the red giant phase lead- 
ing to a C/O ratio in the core that could help to constrain the 
rate of the n C(a, y) 16 reaction. In the present case, the rel- 
evant process is diffusion in the white dwarf envelopes lead- 
ing to a gradually thickening pure He surface layer as the star 
cools. In effect, each of these studies has concentrated on an 
isolated structure in the stellar interior, ignoring the contribu- 
tion of other possible features. Regardless of whether the as- 
teroseismic signal originates from the core or the envelope, 
it is now clear that the relevant structure (C/O profile or He 
surface layer) is very nearly the same in the two stars. The 

bette r absolute qual ity of the fits for GD 358 and CBS 1 14 
I I j 1 I J 

in Metcalfe (2003b) suggest that the C/O profile is the more 
important structure from an asteroseismic standpoint. Indeed, 
this is consistent with the finding by Montg omery et alJ J2003) 
that the pure He surface layer has a smaller effect on the pul- 
sation modes than the uniform He/C envelope, which is itself 
secondary in importance to the C/O profile in the core. 

Of course, there are sound physical reasons to expect com- 
plicated structures in both the core and the envelope of real 
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white dwarfs. We have demonstrated that analyses of two com- 
pletely independent sets of pulsation data with two physically 
distinct classes of models can both lead to reasonably self- 
consistent interpretations of the corresponding internal struc- 
tures. Our challenge is now to synthesize hybrid models that 
can fit the core and envelope structure simultaneously. It is pos- 
sible that such models will exhaust the information content of 
the 1 1 pulsation modes we have observed in these two stars. If 
so, our only recourse is to detect additional modes in GD 358 
and CBS 1 14 and to discover additional DBV stars with rich 
pulsation spectra. Both options will probably require multi-site 
campaigns using 4 m class telescopes since the existing obser- 
vations are already the result of campaigns using 1-2 m tele- 
scopes, and the newly discovered DBV sta rs from the Sloa n 
Digital Sky Survey are significantly fainter (Ni tta et alJ l2005). 
Even so, the future continues to look bright for white dwarf 
asteroseismology. 
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